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Abstract. We present an observational program we started in 1999, to systematically obtain 
mid-resolution spectra of late- type stars, to study in particular chromospheric activity. In par- 
ticular, we found cyclic activity in four dM stars, including Prox-Cen. We directly derived the 
conversion factor that translates the known S index to flux in the Ca II cores, and extend its 
calibration to a wider spectral range. We investigated the relation between the activity mea- 
surements in the calcium and hydrogen lines, and found that the usual correlation observed is 
the product of the dependence of each flux on stellar color, and it is not always preserved when 
simultaneous observations of a particular star are considered. We also used our observations 
to model the chromospheres of stars of different spectral types and activity levels, and found 
that the integrated chromospheric radiative losses, normalized to the surface luminosity, show 
a unique trend for G and K dwarfs when plotted against the S index. 
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1. Introduction 

In 1999 we started a program to systematically obtain spectra of late-type stars, to 
study in particular chromospheric activity. The stars were chosen to cover the spectral 
range from F to M, with different activity levels. Since we were particularly interested in 
the transition to completely convective stars, we included a larger number of M stars in 
our sample. 

Our observations were made at the 2.15 m telescope of the Complejo Astronomico El 
Leoncito (CASLEO), which is located at 2552 m above sea level, in San Juan, Argentina. 
We obtained high-resolution echelle spectra with a REOSC spectrograph. The maximum 
wavelength range of our observations is from 3860 to 6690 A, and the spectral resolution 
ranges from 0.141 to 0.249 A per pixel (R = X/5X w 26400). 

At present, we have about 5500 spectra of 150 stars, ranging from F to M with different 
activity levels. Altough most of the stars are single dwarfs, we have several binaries and a 
few subgiants. Currently, we have four observing runs per year. Details on the reduction 
and calibration procedures can be found in Cincunegui & Mauas (2004). 



2. Flux-Flux calibrations 

Since, unlike most surveys of this kind, the observations in the different spectral fea- 
tures are made simultaneously, our data provides an excellent opportunity to study 



119 



120 



Pablo J.D. Mauas et al. 



-0.5 



-1.5 

















\ 








^# : 

*^ 


- 
- 




\ : 


- - - log 


( C MN) 




log 






— log 




N 


1 







0.5 1 1.5 

B-V 



Figure 1 . Conversion factor between Mount W ilson S index a nd Fhk- The dashed line corre- 



sponds to Noye s et al. (1984) the dotted line t o |Rutten (1984)[ and the full line to our derived 
factor (from Cincuneg ui, Diaz &; Mauas 2007a[ ). 

the correlation between different spectral features and activity indexes. To date, the 
most common indicator of chromospheric activity is the well-known S index, essentially 
the ratio of the flux in the core of the Ca II H and K lines to the continuum nearby 



( Vaughan et al. 1978 ). This index has been defined at the Mount Wilson Observatory, 
were an extensive database of stellar activity has been built over the last four decades. 
However, unlike our program, these observations are mainly concentrated on stars rang- 
ing from F to K, due to the long exposure times needed to observe the Ca II lines in the 
red in faint M stars. For this reason, the S index is poorly characterized for these stars. 
We first obtained the S index for our spectra, integrating with the corresponding profile 
(for details, see Cincunegui, Diaz fc Mauas 2007al). 



The Mount Wilson S index can be converted to the average surface flux in the Ca II 
lines through the relation: 

F H K = ^boil.3410- 4 ^C cf , (2.1) 

where C c f (B — V) is a conversion factor that depends on color. Two different expressions 
are widely used for this factor, the first one given by |Middelkoop (1982)| and corrected 
by Noyes et al. (1984)] and the other one given by Rutten (1984) The deductions used 



in both works to derive C c f involve complex calibration procedures. 

Since we have simultaneous measurements of the S index and the core fluxes, we can 
calculate directly the correction factor as a function of the index and the flux, which 
is shown in Fig. [T] together with the other two expressions (for details and a whole 
discussion, see |Cincuneg ui, Diaz & Ma uas 2007a[ ). 

It is usually accepted that there is a tight relation between the chromospheric fluxes 
emitted in H-a and in the H and K Ca II lines, and that these two features can be 
used to study chromospheric activity. However, most works where this relation has been 
observed found it by using averaged fluxes for both the calcium and the hydrogen lines, 
which were not obtained simultaneously, and were even collected from different sources. 

However, the situation is different when the individual stars are studied separately. In 
Fig. [2] we plot the individual simultaneous measurements of each flux for several stars of 
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Figure 2. H-a vs. Ca II surface fluxes, for st ars of different spectral types, divi ded into 
different color bins, as indicated (from [Cincuneg ui, Diaz Mauas 20Q7a[ ). 

different spectral types and different levels of activity, divided into color bins as stated. 
We also show the linear fits for each star. It can be seen that the behavior is different 
in each case: in some stars both fluxes are well correlated, although the slopes of the 
fits are not the same. In other stars the H-a flux seems to be almost independent of the 
level of activity measured in the Ca II lines, and there are even stars where the fluxes 
are anti-correlated. 

In Diaz, Cincunegui fc Mauas (2007)] we studied the sodium D lines (Dl: 5895.92 A; 



D2: 5889.95 A) in our stellar sample. We found a good correlation between the equivalent 
width of the D lines and the color index (B — V) for all the range of observations. Since 
equivalent width is a characteristic of line profiles that do not require high resolution 
spectra to be measured, this fact could become a useful tool for subsequent studies. 
Finally, we constructed a spectral index (R' D ) as the ratio between the flux in the D lines 
and the bolometric flux. Once corrected for the photospheric contribution, this index 
can be used as a chromospheric activity indicator in stars with a high level of activity. 
Additionally, we found that combining some of our results, we obtained a method to 
calibrate in flux stars of unknown color. 



3. Cycles in M-stars 

One of the main goals of our program was to extend the studies of stellar cycles to M- 
stars, at and beyond the limit for full convectivity. The first star we studied was Proxima 
Centaury, a dMe 5.5 star with strong and frequent flaring activity Cincunegui, Diaz & 
Mauas (2007b). For this star we excluded the spectra taken during flares, computed the 
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Figure 3. Lomb-Scargle periodogram of our observations. The False Alarm Probability levels 
of 50 and 80% are shown, (from |Cincunegui, Diaz fe Mauas 2 007b). 



nightly average of the Ha flux, and calculated the Lomb-Scargle periodogram, which is 
shown in Fig. [3l In the periodogram we found strong evidence of a cyclic activity with 
a period of ~442 days. Similar values for the period were found using three different 
techniques in the time domain (see |Cincunegui, Diaz fe Mauas 2007b| for details). We 
were also able to determine that the activity variations outside of flares amount to 130% 
in S, three times larger than for the Sun. 

Since this star should be fully convective, it cannot support an aQ dynamo, and a differ- 
ent mechanism should be found to explain this result. Recently, |Chabrier fc Kker (2006)] 
showed that these objects can support large-scale magnetic fields by a pure a 2 dynamo 
process. Moreover, these fields can produce the high levels of activity observed in M 
stars (see, for example, IMauas fc Falchi 19961 This a 2 dynamo does not predict a cyclic 
activity. However, our observations suggest that this cool star has a clear period. 

We also studied the spectroscopic binary system Gl 375 (Diaz, Gonzalez, Cincunegui, Mauas 2007). 
We first obtained precise measurements of the orbital period (P = 1.87844 days) and sep- 
aration (a = 5.665 R©), minimum masses and other orbital parameters. We separated 
the composite spectra into those corresponding to each component, which allows us to 
confirm that both components are of spectral type dMe 3.5. 

To study the variability of Gl 375, besides using the spectra obtained at CASLEO, 
we also employed photometric observations provided by the All Sky Automated Survey 
(AS AS Pojmanski 2002). We calculated the Lomb-Scargle periodogram for these data, 
and obtained a a distinct peak corresponding to a period of P p hot = 1.876667 days, 
a period resembling very closely the measured orbital period. We believe this harmonic 
variability is produced by spots and active regions in the stellar surface carried along with 
rotation. This would imply that the rotational and orbital periods are synchronized, as 
is expected for such a close binary. 

To verify if this is indeed the case, we phased the data to the obtained period, for two 
different seasons (Fig. HJ. The sinusoidal shape of the variation is evident, although the 
amplitude of the modulation is different in both cases, probably an indication of different 
area covered by starspots or active regions, and therefore different activity levels, in each 
epoch. Therefore, the amplitude of this modulation can be used as an activity proxy, and 
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Figure 4. AS AS photometry phased t o the orbital period for two different epo chs. Left: 
2002.5-2003.5 Right: 2006. (from |Dfaz, Gonzalez, Cincunegui, Mauas 2007) . 




f K (prim) [ergs cm 2 A 1 s '] 

Figure 5. Comparison of the fluxes in the Ca II K line for both components. The error bars 
correspond to a 10% error in the line flux es, and the dotted line is the identity relation, (from 
|Dfaz, G onzalez 7~Cincunegui, Mauas 2007[ ) . 



indicates that the system exhibits a roughly periodic behavior of 2.2 years (or 800 days). 
The same period was found in the mean magnitude of the system and in the flux of the 
Ca II K line, although the Ca flux variations occur 140 days ahead of the photometric 
ones, a behavior that has been previously observed in other stars (Gr ay et al. 1996[ ). The 
agreement between the behavior of the three observables is remarkable because of the 
different nature of the observations and the different instruments and sites where they 
were obtained. 

Another interesting result of this work is that the activity of Gl 375 A and Gl 375 B, as 
measured in the flux of the Ca II K lines, are in phase, as can be seen in Fig.[5j There is an 
excellent correlation between the levels of chromospheric emission of both components, 
implying a magnetic connection between them. Due to its vicinity and relative brightness, 
this system presents an interesting opportunity to further study this type of interaction. 

We also studied the long-term activity of two other M dwarf stars: Gl 229 A and Gl 
752 A, using again the Ca II - K line-core fluxes measured on our spectra and the ASAS 
photometric data. Using the Lomb-Scargle periodogram, we obtained a possible activity 
cycle of -4 and ~7 yrs for Gl 229 A and Gl 752 A, respectively (jBuccino et al. 201 ip . 
This work was complemented by other studies, were we investigated the presence of 
activity cycles using IUE data ([Buccino fc Mauas 20081 and IBuccino &; Mauas 2 009). 
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Figure 6. Normalized 0mt versus the S index. The squares represent the K star mod- 
els from |Vieytes, Diaz &; Mauas (2009) | and the triangles indicate the G star models from 
Vieytes, Mauas & (Jincunegui (20Q5j[ 

4. Chromospheric modeling 

We also used our spectra as observational basis to model the chromosphere of different 
type of stars ans activity levels, following the path initiated with the model of Ad-Leonis 
in its quiescent state (Mauas & Falchi 1994), which was completed with a study of other 
dM 3.5 stars of different activity levels (|Mauas et al 1997)) . We first computed models of 
the Sun as a star and 9 solar analogues (Vieytes, Mauas & Cincunegui 2005) which were 
followed by models for 6 K stars, "analogs" of epsilon-Eridani (Vieytes, Diaz & Mauas 
2009). For most stars, we built two models, to match the observations in its minimum and 
maximum levels of activity, and found that the differences in the atmospheric structure 
for a star in its maximum and minimum activity levels are similar to the changes seen 
between two different stars. 

We also found that the integrated chromospheric radiative losses <p[ nt , normalized to 
the surface luminosity, show a unique trend for G and K dwarfs when plotted against the 
S index, (see Fig. [6]). This might indicate that the same physical processes are heating 
the stellar chromospheres in both cases. We calculated an empirical relationship between 
S and the energy deposited in the chromosphere, which can be used to estimate the 
energetic requirements of a given star knowing its chromospheric activity level. 
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Discussion 

STEVE SAAR: You mentioned ASAS data in a couple of cases, I would encourage 
you to look at the ASAS data for Proxima Centauri. By an odd coincidence we have 
been looking for an x-ray cycle in that recently and it has what appears to be a good 
photometric period but it is more like 8 years rather than the period you quote here. I 
have been wondering if you see a signal in that range. 

PABLO MAUAS: I was thinking it was time for an update in that star, since we have 
four more years of data and we certainly have to use the ASAS data. We will have a look 
at it. 

JURGEN SCHMITT: Alpha Cen must be perfectly suited for observations you have - 
have you looked at Alpha Cen? 

PABLO MAUAS: Yes, we have, but we had a problem with the gain of the pointing 
screen at the observatory so we had to end those observations three or four years ago. 
But we have a hint of a cycle in the K-star, Alpha Cen B. 

MARK GIAMPAPA: I think it is great that you are getting cycle periods for M dwarf 
stars. Just a comment on the H-a and Calcium K. The models that I have constructed 
and looked at show that they are very segregated in their formation regions and so you 
could have some, at least in models, disjoint behavior, but I am surprised that you don't 
see in observations pretty direct correlation. But, in any event I could imagine in the 
upper atmosphere you could have flare-like behavior that may not propagate into the 
lower atmosphere and that could lead to some deep correlation between Calcium and 
H a. 

PABLO MAUAS: Flare activity we usually can detect because we take two spectra 
to eliminate cosmic rays, so before we do anything with our data we check to see if the 
spectra are more or less the same. So if you had a flare, that wouldn't happen particularly 
not for M stars integration times are very long. So usually we leave out all the flares. 
For Proxima for example that's hard work because you can have four or five flares a day. 
So usually those extremes are eliminated from our data. I think there can be different 
things for these different correlations between Calcium and H-a. One can be filling factor 
problems. Perhaps you can explain the same integrated emission with different filling 
factors and different contrasts but they won't give the same relation between K-a and 
Calcium. The other one is the orientation - it is not the same if you observe the star 
pole-on or from the equator. So we are exploring that, and preliminary results are shown 
in our poster. 
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